Diffusive Spin Dynamics in Ferromagnetic Thin Films with a Rashba Interaction 
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In a ferromagnetic metal layer, the coupled charge and spin diffusion equations are obtained in 
the presence of both Rashba spin-orbit interaction and magnetism. The misalignment between the 
magnetization and the nonequilibrium spin density induced by the Rashba field gives rise to Rashba 
spin torque acting on the ferromagnetic order parameter. In a general form, we find that the Rashba 
torque consists of both in-plane and out-of-plane components, i.e., T = T±y x m + T||rh x {y x rfi). 
Numerical simulations on a two-dimensional nanowire consider the impact of diffusion on the Rashba 
torque and reveal a large enhancement to the ratio T\\/T± for thin wires. Our theory provides 
an explanation for the mechanism driving the magnetization switching in a single ferromagnet as 
observed in the recent experiments. 

PACS numbers: 75.60.Jk 



The manipulation of spin degrees of freedom and spin- 
charge conversion are at the core of the rapid develop- 
ing field of spintronics [1]. A semiconductor-based two- 
dimensional electron gas (2DEG) lacking inversion sym- 
metry is known to electrically generate nonequilibrium 
spin density through the spin-orbit interaction [2 , '3] . The 
same type of spin-orbit interaction, named after Rashba 
[1], is the driving force behind numerous recent interest- 
ing discoveries, e.g., the well-known spin-Hall effect [5|,|6|. 

Besides its dominating role in semiconductor struc- 
tures, Rashba spin-orbit interaction is expected to 
emerge in a thin ferromagnetic layer, such as cobalt (Co), 
sandwiched asymmetrically between a heavy metal thin 
film (Pt) and a layer of metal oxides (AlO^^) [3-0. In 
such a quasi- two-dimensional metal layer (see Fig(T]), the 
effective field Bn generated by the spin-orbit interaction 
is predicted to excite the ferromagnetic order parameter 
by a charge current [ll'-^13^, which has been confirmed 
by several experiments [7-10]. This spin torque, coined 
Rashba torque, falls into to a broader family where the 
spin-orbit interaction enables a transfer of angular mo- 
mentum between the spin and orbital degrees of freedom, 
and has been observed in diluted magnetic semiconduc- 
tors [ii-[3. The same type of spin-orbit-induced torque 
is predicted to improve current-driven domain wall mo- 
tion [12I, [l3| as it has been shown experimentally |iq| . 
Recently, Miron et al., [9] has demonstrated the current- 
induced magnetization switching in a single ferromagnet, 
which represents an outstanding alternative to the cele- 
brated Slonczewski-Berger spin-transfer torque [l3, [HI 
that requires noncohinear magnetic textures such as spin 
valves or domain walls [19]. 

In a nonmagnetic 2DEG, the diffusive spin dynam- 
ics in the presence of Rashba spin-orbit interaction and 
D'yakonov-Perel spin relaxation [20] has attracted sig- 
nificant attention |2ll - l23L resulting in further intrigue 
such as spin-Hall edges |2^. We foresee that in a metal- 
lic ferromagnetic layer accommodating both a Rashba 
spin-orbit interaction and an exchange splitting, the com- 




FIG. 1: A schematic view of a Co film sandwiched by a metal- 
oxide layer (AlO^) and a heavy metal (Pt). The magnetiza- 
tion direction m is arbitrary, the in-plane charge current is 
jfe, and Br is the Rashba effective field perpendicular to je. 



petition between spin relaxation (induced by random 
magnetic impurities and D'yakonov-Perel) and the spin 
precession enforced by both exchange and Rashba field 
gives rise to complex spin dynamics that is important to 
current-driven magnetization manipulation and anoma- 
lous Hall effect. Meanwhile, theoretical accounts on 
Rashba torque are, so far, limited to an infinite medium 
111 ll2|, where diffusive motion is absent. In this Let- 
ter, we employ the Keldysh technique [25] to derive a 
diffusion equation describing the coupled dynamics of 
charge and spin in a quasi-two-dimensional ferromag- 
netic layer submitted to a Rashba spin-orbit interac- 
tion. We show that the coupling between the magnetism 
and the spin-orbit interaction, as mediated by the elec- 
trons through s — d exchange, leads to a Rashba torque 
(T) that has both out-of-plane and in-plane components, 
i.e., T = T±y x rfi -\- T\\rfi x {y x rh). In the case of a 
narrow magnetic wire, we show that the in-plane torque 
T|| can be enhanced significantly. 

The total Hamiltonian for conducting electrons (with 
effective mass m) in a quasi-two-dimensional ferromag- 
netic layer inx — y plane is (fi = 1 is assumed throughout) 



HT = —^a&-{kxz)^ -A^cO- -Tfi^W (1) 
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where k is the momentum, a is the Rashba constant, 
a the Pauh matrices, /S^^c the ferromagnetic exchange 
sphtting, rh the magnetization direction, and W = 
SjLi ~ ^j) the spin-independent impurity poten- 
tial with centers located at Rj. An electric field is ap- 
plied and to be included in the following discourse. Using 
Dyson equation [25], the advanced and retarded Green's 
functions and G^, Keldysh function (5^, and the 
self-energ y f ]"^^^^^ are related by the quantum kinetic 
equation ^ ^ 

^qR^-Iq^ _ = t^G^ _ G^s^, (2) 

where all Green's functions are the full functions. 

To obtain a diffusion equation from Eq.(j2j), we as- 
sume short-range (5-function type impurity scatterers at 
low concentration with a weak coupling to electrons 
[23| . then a second-order Born approximation is justified, 
i.e., I]^'^'^(r,rO = 6{r,r')G^^^^^ {r,r)/{mr). The 
momentum relaxation rate 1/r due to spin-independent 
impurities is evaluated at the Fermi energy. The qua- 
siclassical distribution function g = Qk^^iT^R)^ as the 
Wigner transform of the Keldysh function {r^ t; r', t')^ 
is obtained by integrating out the relative spatial- 
temporal coordinates while retaining the center-of-mass 
ones i^ = (r + r')/2 and T = (t + t')/2. A gradient 
expansion on both sides of Eq.(|2]) followed by a Fourier 
transform in time domain gives us 
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where p is the density matrix. In Eq.(|3]), Q = u -\- i/r, 
bk = \Axc'n^/2-{-ak x the operator U = & -{A^ 01^^/2 -\- 
ak X z)/bk, and {•, •} denotes the anticommutator. 

We solve Eq.(|3]) according to the discussion outlined 
in Ref.[23^: find a formal solution g = ^[K^^^ + K^^^] in 
terms of the lowest order approximation 



G^{k,e)p{e)-p{e)G^{k,e) 



and a higher-order gradient correction 
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that is accounted by perturbation. We then apply the 
formal solution in K^^"^ to obtain the zeroth order ap- 
proximation that is to be substituted into K^^"^ to obtain 
a gradient correction, thus arriving at the first order ap- 
proximation to a. W e repeat this procedure to the sec- 
ond order for g [2^. This gradient expansion scheme is 
appHcable as long as the spatial gradient of the quasi- 
classical distribution function is smooth at the scale of 
Fermi wavelength, i.e. dr <^ kp. 



After angle averaging in momentum space and Fourier 
transforming the second-order approximation back to 
real time, we have a diffusion- type equation for the den- 
sity matrix 

d . 1 . 1 1 

— p + — p - - — {z X a) ■ p{z X 0-) + :^{crmp(Tm - P) 

= DV'^p + iC[z X 6-, Vp] - B{z X (T, Vp} - zA^c[<3-m, p] 
+ r[(rh X V)zp - Vp - {z x a) - {z x cr) - V p&rn] 
- 2R{a^, (rh x V),p} (6) 

where am = & ■ rfi and we concentrate on quantities at 
Fermi energy. The diffusion constant is D = rVp/2^ given 
vf the Fermi velocity. A^c = (Aa;c/2)/(4^^ + 1) where 
= (A^^/4 + a^A:|.)r^. The other parameters in Eq.® 
are 
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where we identify r^^c as the D'yakonov-Perel relaxation 
time. Equation (jG]) is valid in the dirty limit (i.e. ^ 1), 
which permits the approximation 1 + 4^^ ~ \ throughout 
the following discussion. We decompose the density ma- 
trix as /3 = n/2 + S ■ <T to introduce the charge density n 
and the spin density S. Spin transport in ferromagnetic 
layers in a real experimental setup [7, lo] is exposed to 
random magnetic scatterers, for which an isotropic spin- 
flip relaxation S/Tsf is introduced phenomenologically. 
In total, we have 

at 

+ rV^ • rhn + RVz • rh{S • rh), (7) 

at T\\ " 

- AxcS X rh — 7p— m x (S' x rh) 

xc 

-h BV^n + 2CVz X 2R{rh • V^n)rh 

+ F [rh X ( X S') + X (rh X S)] , (8) 

where = £ x V. Accounting for both D'yakonov- 
Perel and random magnetic impurities, rate l/ry = 
^/txc + ^/^sf measures the relaxation of the spin den- 
sity S'li = SxX + Syy^ as l/r^ = 2/txc ^ ^/Tsf does to 
= SzZ. Parameter Txc sets a time scale for the decay 
of the transverse (to rh) component of the spin density, 
thus contributing directly to the spin torque [26]. 

For a broad range of the ratio akp / A^c^ Eqs.([7j) and 
([8]) govern the full spin dynamics in a ferromagnetic 
layer and constitute the main result of this Letter. We 
can readily show [13] that the absence of magnetism 
{Axe = 0) leaves the B term a source generating spin 
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density electrically [3|, |23| . When spin-orbit coupling van- 
ishes {a = 0), the first two lines in Eq.(j8j) survive to 
describe a diffusive motion of spin density in a ferromag- 
netic metal, which agrees excellently with early results in 
the corresponding limit jl^, HsJ • The C term represents 
a coherent precession of the spin density around the ef- 
fective Rashba field. The precession of the spin density 
(induced by the Rashba field) around the exchange field is 
described by the F term, is thus at a higher order (com- 
pared to C) in the dirty limit for F = AxctC/2. The 
R-teim contributes to a magnetization renormalization. 

To illustrate the spin dynamics and the Rashba spin 
torque embedded in Eq.(j8]), we restore the electric field 
of strength (o applied along the x direction by a shift 
V ^ V + eS'xde [iH. This Letter is focusing on a 



weak spin-orbit coupling akr < ^xc [Ul, ll2|; i-e., the 
spin density is aligned dominantly along the local mag- 
netization. Thus, the deviation due to effective Rashba 
field Br (X je X z (along the y direction in the present 
setup) is considered as a perturbation. We may ap- 
proximate the energy derivative d^S ^ PpMp'rh and 
deU ~ np/^F = A/f, given Mf the density of states and 
Pf the polarization; both quantities are at Fermi energy 

Consider an infinite homogeneous ferromagnetic layer 
we replace V by 6(^x8^ and Eq. ([8]) is 



— - = S m X [S X m) S x m 

Ot Tsf Txc TA 

+ e^A^F [^PfCv xrfi^ PFTrfi x {y x rfi)] . (9) 

where the higher-order contribution proportional to R is 
discarded and ta = 1/ ^xc characterizes the time scale of 
the precession of the spin density around the magnetiza- 
tion. The spin torque exerted on the local magnetization 
by the nonequilibrium spin density is given by 



— S X rh 
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taking into account a fieldlike spin precession and de- 
phasing of the transverse component that is essentially 
the Slonczewski-Berger type spin-transfer torque. At a 
stationary state dS/dt = 0, Eq.(j9]) leads to a general 
form T = T±y x rh -\- T\\rfi x {y x rh) that is 
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[{2r]C + PT)y x rh 

+(7?F - 2/3C)m X {yxrh)], (11) 



where je = e^UFrS /m is the current density, C = 
ta(1/ts/ + 1/T^c), = 1 + C^a/Txc, and p = ta/tsJ 
[28|. The first term in Eq. (pT]) is an out- of -plane torque 
driven by the effective field Br and is perpendicular to 
the (m, Br) plane. Using the data from experiments 
0, we take a ^ 10~^^ eVm, polarization Pp = 0.5, sat- 
uration magnetization Mg ~ 1.09 x 10^ A m~^, and the 



current density je = 10^ A cm~^, the effective field gen- 
erating the out- of -plane torque is estimated to be of the 
order of 1 T. 

The second in-plane torque arises from the change of 
the spin density induced by its precession around the 
exchange field (ex F) as well as to the presence of spin flip 
in the layer (ex /3), inducing a spin density that is parallel 
to Br X rh {ot y X rh) [2?], hence leading to a torque 
in the (rh, y) plane. In other words, the in-plane torque 
is driven by an effective field along the direction Br x 
rh. In the limit Tsf oo, the ratio T\\/T±_ ^ 
approaches a constant. When considering T^c oo and 
F <C C (dirty limit), the in-plane torque diminishes and 
Eq. (fTT]) is in a good agreement with Ref. fn\. 

The impact of diffusion on the Rashba spin torque is 
appreciated by numerically solving Eq. (j8j) in a semi- 
infinite ferromagnetic wire with a width L along the y 
axis. The wire is open on the transport direction x, along 
which an electric field is applied (thus V^^ eS'xde)- To 
compare with the recent experiment 0, the magnetiza- 
tion is aligned perpendicularly along z. Two boundary 
conditions are available in the transversal direction {y): 
vanishing spin density S(ii = ^L/2) = or vanishing 
spin current Jy{y = i^/2) = at the boundaries. The 
former condition indicates an absorbent surface where 
spin-flip relaxation at the boundaries is fairly strong, thus 
suppressing the spin density. A vanishing spin current 
designates a fully reflective boundary where the spin cur- 
rent is reflected. In the present simulation, we choose to 
adopt the boundary condition according to the work of 
Kato et al., j29|: optical detection on the spatial profile 
of the spin density induced by the spin-Hall effect in a 
n-GaAs 2DEG reveals that the spin density vanishes at 
the interfaces. We believe that the recent experiments 
on Rashba spin torque in ferromagnetic layers [^3, [sl, [lo| 
fall into this picture. 

Figure [2] shows the spatial profile of the three compo- 
nents of the spin density: Sx (panels a and d), Sy (panels 
b and e), and Sz (panels c and f), for various magnitudes 
of a. The spin-orbit interaction drives the spin dynamics 
in two prominent ways: the first is explained by Edelstein 
[3]; i.e., a nonequilibrium spin density at y direction is 
generated electrically, becoming a source to a nonvanish- 
ing Sy in the center of the wire. Second, the spin density 
precesses around the total field combining the exchange 
{/S.xc'rh X S) and the effective fields {2CV z x S). In the 
case of a wide magnetic wire, this precession is local- 
ized near the edges of the wire and vanishes towards the 
center where the nonequilibrium spin density generated 
by spin-orbit coupling gets robust along the Rashba field 
{y direction). Consequently, the Sx component, though 
much smaller in magnitude than Sy, is peaked at the 
boundaries; see Figsl2](b) and[2fc). The profile of the Sz 
component has an opposite sign at two edges, which is 
driven by the rotation of the spin accumulation around 
the effective Rashba fields (pointing oppositely to each 
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FIG. 2: (Color online) Spatial profile of nonequilibrium spin 
density Sx (a),(d), Sy (b),(e), and Sz (c),(f) for a width L = 
50nm (a),(b),(c) and L = lOnm (d),(e),(f) for a = 0.1 - 
1 X 10"2 eV •nm. The parameters are r = 10 ^^s, exchange 
splitting Axe = lO^^s"^, spin relaxation time Ts/ = 10~^^s, 
diffusion constant D = 125 x 10~^m^ s~^. 



other) generated while diffusing transversally along op- 
posite directions (i.e., negative y and positive y axes). 

In the case of a weak spin-orbit coupling addressed in 
this Letter, the spin density profile is mainly symmetric 
along the (transverse) y direction. In the opposite limit of 
a strong spin-orbit couphng, the spin-Hall effect induces 
a large spin imbalance along the transverse direction of 
the wire, which favors a deeper asymmetry in the spatial 
profile (not shown here). The narrow width (shorter than 
the spin- flip diffusion length defined by r^/) of the wire 
prevents the relaxation of the spin density such that its 
effective magnitude in the middle of the wire may differ 
from its value in an infinite medium. 

Therefore, one expects a strong influence of the wire 
width on the magnitude and sign of the Rashba spin 
torque. Figures [Sfa) and [3jb) display the magnitude of 

L /2 

the average torque J_j^^^Tdy /L at different wire widths 
as a function of a. For a fixed a, by decreasing the wire 
width the in-plane torque Ty changes its sign and the 
magnitude can be enhanced dramatically, while the mag- 
nitude of the out-of-plane torque Tj_ is just moderately 
reduced. Figure [3l^c) shows the ratio T\\/T± as a function 
of the wire width for different spin diffusion lengths. Sev- 
eral features deserve attention, (i) For a given spin-flip 
relaxation time (therefore the spin diffusion length), the 
ratio approaches a constant when increasing the width of 
the wire, while its absolute value increases exponentially 
when decreasing the width, (ii) For a wide wire with a 
large spin diffusion length, the ratio tends towards the 
bulk value T\\/T± = Axct/4:. Using the present value 
Axe = 10^^ s~^ and r = 10~^^ s, the bulk value is about 
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FIG. 3: (Color online) Out-of-plane (a) and in-plane (b) 
torques as a function of a for different wire widths, (c) Ratio 
T\\/T± as a function of the wire width for different spin dif- 
fusion lengths. The other parameters are the same as in Fig. 

n 



0.03, which agrees well with the numeric value (black 
solid line) in Figl3] (c). (iii) The sign of the ratio can be 
switched by the wire width or the spin relaxation rate, 
as is reflected in the second term in Eq. (pT]) . 

In a recent work, Miron et al, [9] has demonstrated a 
current-induced magnetization reversal in a perpendicu- 
larly magnetized single ferromagnet. In Eq. (pT]) of this 
Letter, the effective field (producing the in-plane torque) 
parallel to Br x rh explains the effective perpendicular 
magnetic field Bs^ required to accomphsh the magne- 
tization switching [9]. In the Co layer as sandwiched 
by metal oxides and heavy metal with a large spin-orbit 
coupling, the spin-flip relaxation length is believed to be 
much shorter than its value in the bulk Co, hence further 
enhancing the in-plane torque: see Figl3](c). If we take 
A^c ^ 5 X 10^"^ s-^ (or 0.3 eV) [30], r,/ ^ 10"^^ s, and 
r ^ 10-^^ s, then we have p ^ 0.002, C ~ ta/T^c ~ 0.5, 
and 77 is of order 1. A quick estimate suggests that 
the perpendicular switching field generating the in-plane 
torque is about 100 mT, which agrees with the estima- 
tion in Ref.[9]. In the case of magnetic domain walls, 
the present torque acts like a transverse field that can 
increase the Walker breakdown limit and enhance the 
range for the nonadiabatic current-driven domain wall 
motion [31]. Miron et al., exploited this characteristic to 
interpret large current-driven domain wall velocities in 
perpendicularly magnetized domain walls [lo| . 

Most of the previous works on transitions metals have 
been carried out in fairly complex structures where inter- 
facial Rashba spin-orbit coupling has not been evaluated 
experimentally or theoretically. Identifying the relevant 
mechanisms is an urgent need [32]. Among them, the 
microscopic description of interfacial Rashba spin-orbit 
coupling in ultrathin layers as well as the spin dynamics 



5 



in the presence spin-orbit coupling will provide essential 
information to the electrical manipulation of spins in low 
dimensional systems. 
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